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ABSTRACT 

Cosmic rays produced in cluster accretion and merger shocks provide pressure to the 
intracluster medium (ICM) and affect the mass estimates of galaxy clusters. Although 
direct evidence for cosmic-ray ions in the ICM is still lacking, they produce 7-ray 
emission through the decay of neutral pions produced in their collisions with ICM 
nucleons. We investigate the capability of the Gamma-ray Large Area Space Telescope 
(GLAST) and imaging atmospheric Cerenkov telescopes (lACTs) for constraining the 
cosmic-ray pressure contribution to the ICM. We show that GLAST can be used to 
place stringent upper limits, a few per cent for individual nearby rich clusters, on 
the ratio of pressures of the cosmic rays and thermal gas. We further show that it is 
possible to place tight (< 10%) constraints for distant {z < 0.25) clusters in the case 
of hard spectrum, by stacking signals from samples of known clusters. The GLAST 
limits could be made more precise with the constraint on the cosmic-ray spectrum 
potentially provided by lACTs. Future 7-ray observations of clusters can constrain the 
evolution of cosmic-ray energy density, which would have important implications for 
cosmological tests with upcoming X-ray and Sunyacv-Zcl'dovich effect cluster surveys. 
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— radiation mechanisms: non-thermal — gamma-rays: theory. 
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1 INTRODUCTION 

Clusters of galaxies are potentially powerful observational 
probes of dark energy, the largest energy budget in the Uni- 
verse causing the cosmic acceleration (e.g., Haiman, Mohr, & 
Holder 2001; Albrecht et al. 2006). Most of the cosmological 
applications using clusters rely on the estimates of their to- 
tal virial mass — quantity which is difficult to measure accu- 
rately in observations. Clusters offer a rich variety of observ- 
able properties, such as X-ray luminosity and temperature 
(e.g., Rosati, Borgani, & Norman 2002), Sunyaev-Zel'dovich 
effect (SZE) flux (e.g., Carlstrom, Holder, & Reese 2002), 
gravitational lensing of distant background galaxies (e.g.. 
Smith et al. 2005; Dahle 2006; Bradac et al. 2006), and ve- 
locity dispersion of cluster galaxies (e.g., Becker et al. 2007) 
and proxies for cluster mass. 

One of the most widely used methods for measuring 
cluster masses relies on the assumption of hydrostatic equi- 
librium between gravitational forces and thermal pressure 
gradients in the intracluster medium (ICM) (Sarazin 1986; 
Evrard, Metzler, & Navarro 1996). Current X-ray and SZE 
observations can yield mass of individual clusters very pre- 
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cisely based on accurate measurements of the density and 
temperature profiles (Pointecouteau, Arnaud, & Pratt 2005; 
Vikhlinin et al. 2006; LaRoque et al. 2006). However, the ac- 
curacy of the hydrostatic mass estimates is currently limited 
by nonthermal pressure provided by cosmic rays, turbulence, 
and magnetic field in the ICM (Ensslin et al. 1997; Rasia 
et al. 2006; Nagai, Vikhlinin, & Kravtsov 2007a, and ref- 
erences therein). This nonthermal bias must be understood 
and quantified before the requisite mass measurement accu- 
racy is achieved. Comparisons with the mass estimates from 
gravitational lensing can provide potentially useful limits on 
this nonthermal bias (see e.g., Mahdavi et al. 2007). How- 
ever, present observations do not yet constrain the nonther- 
mal pressure in the regime in which it dramatically affects 
the calibration of the hydrostatic mass estimates. If not ac- 
counted for, these nonthermal biases limit the effectiveness 
of upcoming X-ray and SZE cluster surveys to accurately 
measure the expansion history of the Universe. Detailed in- 
vestigations of sources of nonthermal pressure in clusters are 
thus critical for using clusters of galaxies as precision cos- 
mological probes. 

There is growing observational evidence for the non- 
thermal activity in clusters. For example, radio and hard X- 
ray observations of clusters suggest presence of relativistic 
electrons. This also implies presence of relativistic protons 
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produced in the shock that accelerated these electrons. How- 
ever, the signature 7-ray emission due to decays of neutral 
pions produced in the collisions of cosmic rays with nucle- 
ons in the ICM has not been detected. From non-dctcctiori 
of 7-ray emission from clusters with the Energetic Gamma 
Ray Experimental Telescope (EGRET) in the GeV band 
(Reimer et al. 2003; but see also Kawasaki & Totani 2002; 
Scharf & Mukherjee 2002), constraints have been placed on 
the fraction of cosmic-ray pressure in nearby rich clusters at 
less than ~20% (Ensslin et al. 1997; Miniati 2003, Virgo and 
Perseus clusters) and less than ~30% (Pfrommer & Enfilin 
2004, Coma cluster). Similar constraints arc obtained using 
the Whipple Ccrcnkov telescope in the ToV band (Perkins 
et al. 2006). These measurements indicate that the cosmic 
rays provide relatively minor contribution to the dynamical 
support in the ICM (e.g., Blasi 1999). However, the current 
constraints are too loose for the future cluster-based cosmo- 
logical tests. 

The next generation of 7-ray detectors, such as Gamma- 
ray Large Area Space Telescope (GLAST) and imaging at- 
mospheric Cerenkov telescopes (lACTs), will be able to pro- 
vide dramatically improved constraints on the cosmic-ray 
pressure in clusters, and may even detect 7-ray radiation 
from several rich clusters (Ando et al. 2007a, and references 
therein). The Gij45r satellite, which is soon to be launched, 
is equipped with the Large Area Telescope (LAT) that en- 
ables all sky survey with GoV 7-rays. Several lACTs arc cur- 
rently working or planned for detecting TeV 7-rays, which 
include HESS, MAGIC, VERITAS, and CANGAROO-III. 
Confronting the recent advances in 7-ray astronomy as well 
as growing interests in dark energy studies, in the present 
paper, we investigate the sensitivity of these detectors to 
high-energy 7-rays of cosmic-ray origin. 

We first show updated sensitivities of GLAST and 
lACTs for nearby rich clusters following Pfrommer & Enfilin 
(2004). In particular, GiAST would be able to constrain the 
cosmic-ray energy density in such clusters to better than a 
few per cent of the thermal energy density, while lACTs 
would be useful to constrain the cosmic-ray spectrum. We 
then consider stacking many 7-ray images of distant clusters 
to probe the evolution of cosmic-ray pressure. We show that, 
by stacking many massive clusters, the upcoming GLAST 
measurements will have the statistical power to constrain 
the cosmic-ray pressure component to better than ~10% of 
the thermal component for clusters out to z < 0.25. These 
forthcoming measurements will be able to place stringent 
limits on the bias in the cluster mass estimates and hence 
provide important handle on systematic uncertainties in cos- 
mological constraints from upcoming X-ray and SZE cluster 
surveys. 

Throughout this paper, we adopt the concordance cos- 
mological model with cold dark matter and dark energy 
(ACDM), and use Q.^ = 0.3, Qa = 0.7, Hq = 100 h km 
s"^ Mpc"^ with h = 0.7, and as = 0.8. 



2 7-RAY PRODUCTION DUE TO 
PROTON-PROTON COLLISIONS 

Cosmic-ray protons are injected in the ICM through the 
shock wave acceleration, and the momentum distribution 
follows the power law, pp with ap ~ 2-3. These cosmic- 



ray protons then interact with the surrounding ICM (mostly 
nonrelativistic protons), producing neutral and charged pi- 
ons; the former decays into two photons (tt" — > 27) while the 
latter into electrons, positrons, and neutrinos. The volume 
emissivity of the 7r°-decay 7-rays (number per volume per 
unit energy range) at distance r from the cluster center is 
given as (e.g., Blasi & Colafrancesco 1999) 



poo poo 

2nH(r)c / dE^ / dE, 



dap 



dE^ 



{En, Ep) 



Ep,th(E^) 

np{Ep,r) 
VE-i-mic^' 



(1) 



where mn and E^ is the mass and energy of the neutral pion, 
-E,r,min = E + m% / 4,E IS the minimum pion energy required 
to produce a photon of energy E, and similarly Ep^th{En) 
is the minimum energy of protons for pion production. The 
density of ICM, nYi{r), is very well measured by the X-ray 
observations of bremsstrahlung radiation from thermal elec- 
trons, and the cross section of the proton-proton collision 
for pion production, dapp/dE-n-, can be calibrated using lab- 
oratory data. The distribution function of cosmic-ray pro- 
tons np(Sp,r) depends on the injection power, spectrum, 
and spatial distribution of cosmic rays. By specifying these 
ingredients, we can predict the 7-ray flux from a cluster. 

In practice, we use a fitting formula as well as cluster 
parameters given in Pfrommer & EnBlin (2004); for the for- 
mer, we briefly summarize it in Appendix A. In addition, one 
should also note that electrons and positrons produced by 
TT* decays can scatter cosmic microwave background (CMB) 
photons up to 7-ray energies. For a while, we neglect this 
secondary process, but revisit it in Section 6 and show that 
it is in fact negligible under most of realistic situations. 



2.1 Cosmic-ray power and spectrum 

The cosmic-ray pressure Pp and energy density pp, which are 
the quantities that we want to constrain, are directly related 
to the injection power of the cosmic rays. The cosmic-ray 
spectrum is measured to be a power law with the index of 
ap = 2.7 in our Galaxy, but in the clusters it is perhaps 
harder, since they can confine cosmic rays for cosmological 
times (Volk, Aharonian, & Breitschwerdt 1996; Ensslin et 
al. 1997; Bcrczinsky, Blasi, & Ptuskin 1997). We thus adopt 
harder spectrum with ap = 2.1 and 2.4, but also use ap = 
2.7 as a limiting case. 

It is also possible that the injection of the cosmic rays 
and thus their energy density pp arc intermittent. Although 
it is interesting to constrain the source property by measur- 
ing such 7-ray variability, this is not the primary focus in 
the present paper. Instead, we concentrate on constraining 
energy density pp averaged over GLAST exposure time. For 
the sensitivity of GLAST, we consider the result of one- year 
all-sky survey, which corresponds to ~70-day exposure to 
each source as the field of view is ~20% of the whole sky. 
Therefore, any time variability within this 70-day duration 
is smeared out. 
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2.2 Radial distribution 

We define quantities Xp and Yp as ratios of energy density 
and pressure of cosmic rays to those of thermal gas, respec- 
tively, i.e., 

Xp^-^, Y,^^. (2) 
Pth ^th 

In general, these depend on the radius, but the concrete de- 
pendence is totally unknown. Various mechanisms supplying 
the cosmic-ray protons have boon proposed, which produce 
characteristic and diverse profiles of Xp and Yp. We thus 
parameterize them using a simple power-law: 

Xp{r) = XpiRsoo) (^^y , 

Yp{r) = Yp{R,oo) (^-^y , (3) 

where i?A (here A = 500) is the radius at which the en- 
closed spherical overdensity is A times the critical density 
of the Universe at the cluster's redshift,^ where the cluster 
mass A/a is traditionally defined with the X-ray and SZE 
measurements. We note that this approach ignores boosts 
in 7-ray flux caused by dumpiness. The constraints derived 
using a smooth model hence provide a conservative upper 
limit on Xp and Yp. 

Wc first focus on Xp, and later discuss Yp. The rela- 
tion between 7-ray intensity and Xp is summarized in Ap- 
pendix A and that between Xp and Yp is discussed in Sec- 
tion 5. We shall study the dependence of results on (3, for 
which we adopt 1, 0, and —0.5. Below, we outline several 
models that motivate these values of /J. 

2.2.1 Isobaric model 

The simplest model is based on the assumption of /3 = 0, 
i.e., the energy density of cosmic rays precisely traces that of 
thermal gas everywhere in the cluster. The latter is propor- 
tional to temperature times number density of the thermal 
gas, both of which are very well measured with X-rays for 
various nearby clusters. The gas density profile is nearly con- 
stant within a characteristic core radius rc, beyond which it 
decreases as a power law, while temperature profile is almost 
constant. The core radius and outer profile are rc = 300 kpc, 
r"^-^ (Coma), rc = 200 kpc, r"^ '^ (Perseus), and rc = 20 
kpc, r"^" (Virgo) (see Table 1 of Pfrommer & Enfihn 2004, 
for a more comprehensive list). The latter two clusters have 
an even smaller 'cool core', but this structure gives only a 
minor effect on the 7-ray fiux. 

2.2.2 Large-scale structure (LSS) shocks 

The formation of galaxy clusters is due to merging or accre- 
tion of smaller objects. When this occurs, the shock waves 
are generated at the outskirts of the clusters, somewhere 
around ~3 Mpc from the center, where protons and electrons 

1 We use RAhE(z) = rziT^pec/^ keV)''/^, where rs = 
0.792 /i-i Mpc and ?? = 1.58 for A = 500, rs = 0.351 h"! Mpc 
and ri = 1.64 for A = 2500, and E'^(z) = ^^{1 + 2)* + Ha for 
the flat ACDM cosmology (Vikhlinln et al. 2006). 



are accelerated to relativistic energies (e.g., Loeb & Waxman 
2000; Miniati 2002; Keshet et al. 2003; Gabici & Blasi 2003). 
Unlike electrons that immediately lose energies through syn- 
chrotron radiation and invcrse-Compton (IC) scattering off 
CMB photons, protons are hard to lose energies, and they 
are transported efficiently into the cluster center following 
the motion of ICM gas (Miniati et al. 2001). In order to pre- 
dict the eventual profile of the cosmic-ray energy density, 
one needs to resort to numerical simulations. The recent 
radiative simulations by Pfrommer et al. (2007) show some- 
what jagged shape for the Xp(r) profile, which implies large 
clumping factor. Here, wo model its global structure with a 
smooth profile with /3 — —0.5, ignoring the effects of dumpi- 
ness. On the other hand, they also performed nonradiative 
simulations which rather imply (3=1 profile. Although the 
latter may not be realistic, the effects of cooling and heating 
in clusters arc also somewhat uncertain. Thus, we still adopt 
this model, treating it as an extreme case. 

2.2.3 Central point source 

A central powerful source such as active galactic nuclei or 

cD galaxy might be an efficient supplier of the cosmic rays, 
which diffuse out from the central region after injection. The 
profile of cosmic-ray energy density is r^^ , but truncated at 
a radius that is fax smaller than R^oo for relevant energies 
(Berezinsky et al. 1997; Colafrancesco & Blasi 1998). The 
actual 7-ray detection might therefore cause significant over- 
estimate of the cosmic-ray pressure; we address this issue in 
Section 3.2. 

Numerical simulations of jets from active galactic nuclei 
suggest that temporal intermittency and spatial structure 
might be complicated (e.g., O'Neill et al. 2005). Neither of 
these, however, affect our results that depend on global and 
time-averaged properties. 



3 COSMIC-RAY ENERGY DENSITY IN 
NEARBY GALAXY CLUSTERS 

3.1 Constraints from entire region of clusters 

We first discuss the case of the Coma duster, focusing on 
the region within -R500 = 2.1 Mpc and assuming the iso- 
baric distribution of the cosmic-ray energy density (/? = 0) . 
Fig. 1 shows the integrated 7-ray fiux with photon energies 
above i?min, F{> Smin), for Xp = 0.1. This flux is to be 
compared with the sensitivities of GLAST and lACTs, for 
which one has to take the source extension into account. 
Indeed, the radial extension of the Coma cluster R500 corre- 
sponds to S500 = 1^2, which at high energies exceeds the size 
of the point spread function (PSF), SOpsf{E). We obtain the 
fiux sensitivity for an extended source from that for a point 
source by multiplying a factor of max[l, S5oo/(56psF(i5min)], 
if the sensitivity is limited by backgrounds. On the other 
hand, if the expected background count from the cluster re- 
gion is smaller than one, which is the case for GLAST above 
~30 GeV, the sensitivities for a point source and an ex- 
tended source are identical. The region ~2-30 GeV is where 
the expected background count is smaller than one from the 
PSF area but larger than one from the entire cluster. We as- 
sume that lACTs are limited by background over the entire 



© 2008 RAS, MNRAS 000, 1-12 



4 S. Ando and D. Nagai 



Table 1. Sensitivity to Xp(i?5oo) of GLAST (i?min = 100 MeV) for various values of spectral index of cosmic rays ap, and isobaric and 
LSS shock models for radial distribution or /3. The limits on Xp are set by the 7-ray flux from a region within whichever of the larger 
between the point spread function (59psF(Emin) ~ 3° and the source extension 6500 ■ 
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Figure 1. Flux of 7-ray emission from the region within R500 = 
2.1 Mpc of the Coma cluster, for the isobaric model with Xp = 
0.1 (labeled as 'vr^-decay'). The spectral index of the cosmic-ray 
protons is ap = 2.1 (solid), 2.4 (dashed), and 2.7 (dot-dashed). 
The sensitivity curves of GLAST and lACTs are for a source 
extended by 6500 = 1? 2 (corresponding to if 500 ) 1 while the point- 
source sensitivity of GLASTis also shown as a short dashed curve. 
Flux due to IC scattering and nonthermal bremsstrahlung is also 
shown (dotted; from Reimer et al. 2004). 



energy region, and we multiply the point source sensitiv- 
ity by 6^00 / 56ps¥ with (50psf = 0?1; this is consistent with 
Aharonian et al. (1997) for relevant energy regime. A more 
detailed derivation of this sensitivity is given in Appendix C. 

We also show flux of IC scattering and bremsstrahlung 
radiations from electrons primarily accelerated in the shocks 
(Reimer et al. 2004). The authors suggested that these 
electron components would always be below the GLAST 
and lACT sensitivities, based on constraints from radio, 
extreme-ultraviolet (EUV), and hard X-ray observations. If 
this is the case, the 7-ray detection would imply existence 
of cosmic-ray protons, and be used to constrain the pressure 
from this component (see also, Enfilin, Lieu, & Biermann 
1999; Atoyan & Volk 2001). We give more detailed discus- 
sions about IC mechanisms in Section 6. 

Fig. 1 shows that 7-rays from tt" decays are detectable 



for Xp = 0.1. In particular, the models with different val- 
ues of Dtp predict similar amount of 7-ray fluxes for low- 
energy thresholds {E-mi-n < 1 GeV); GLAST measurements 
can therefore provide constraints on Xp, almost independent 
of ap. Non-detection with GLAST from these nearby clus- 
ters is also very interesting as it provides very tight upper 
limit to the cosmic-ray energy density in clusters. The fluxes 
above ~1 TeV, on the other hand, depends very sensitively 
on Op; lACTs will thus constrain the spectral index. 

In Table 1, we summarize the sensitivity to Xp(i?5oo) 
for GLAST in the case of Emin = 100 MeV, for several values 
of Op and different models of radial distribution of cosmic- 
ray energy density. We also performed the same analysis for 
other nearby rich clusters (Perseus, Virgo, Ophiuchus, and 
Abell 2319), and report their results as well. This indeed con- 
firms that the GLAST constraints on Xp depend only weakly 
on the assumed spectral index. ^ The constraints improve for 
smaller values of /3. For /3 5S 0, the GLAST non-observation 
can place tight upper limits on the cosmic-ray energy den- 
sity at a few per cent level. Even in the case of nonradiative 
LSS shock model (/3 = 1) the constraint is still as good as 
~10% for the Coma. This is a dramatic improvement from 
the EGRET bounds (see, e.g.. Table 3 of Pfrommer & Enfilin 
2004), by more than an order of magnitude. 

On the other hand, the lACT constraints on Xp (with 
Emin ~ 1 TeV) for the Coma cluster and /3 = profile 
are 0.37, 2.3, and 42 for ap = 2.1, 2.4, and 2.7, respec- 
tively. Thus, lACTs will therefore provide constraints on 
the spectral index, which is directly related to astrophysi- 
cal mechanisms of particle acceleration. A similar trend can 
be found in Table 6 of Pfrommer & Enfilin (2004); however, 
the authors apphed point-source flux limit to the (extended) 
clusters and obtained much more stringent sensitivities than 
ours. 

3.2 Direct constraint from large radii 

So far, we treated all clusters but Virgo as point sources. 
Although we showed that the dependence on the assumed 
radial proflle was reasonably weak, a more general approach 
would be to use the resolved image. This is particularly use- 
ful, if the radial proflle cannot be simply parameterized (see 
Section 2.2). Because we are interested in the cosmic-ray 

^ Note that the sensitivity peaks at ap = 2.5. This is because for 
even larger Op, the contribution from low- momentum protons to 
the energy density becomes more significant, while they do not 
produce 7-rays. 



© 2008 RAS, MNRAS 000, 1-12 



Cosmic-ray pressure in galaxy clusters 5 



pressure at -R500 and the 7-ray yields would rapidly drop 
with radius, we here consider constraints in a projected ra- 
dial shell between 612500 and 6300- We mainly focus on the 
Perseus, and assume ap — 2.1; in this case, ^2500 ~ 0?65. 
In order to resolve the inner region, we consider the en- 
ergy threshold of 0.6 GeV, above which the GLAST reso- 
lution becomes smaller than ^2500- The GLAST flux lim- 
its for the outer region and for E > 0.6 GeV correspond 
to the following limits on the fractional energy density: 
Xp,iim{R50o) = 0.099, 0.089, and 0.080, for /3 = 1, 0, and 
—0.5, respectively, which are still reasonably small. In ad- 
dition, these are much less sensitive to the assumed profile, 
thus applicable to more general cases including the central 
source model. The similar procedure predicts sensitivities for 
other clusters: Xp(i?50o) = 0.42 (Coma), 0.14 (Virgo), 0.41 
(Ophiuchus), and 0.55 (Abell 2319), in the case of /3 = 
and Qp = 2.1. Although it is limited to nearby clusters, such 
analysis provides an important handle on the radial distri- 
bution of cosmic-ray ions in clusters. 



4 EVOLUTION OF COSMIC-RAY ENERGY 
DENSITY 

While we could obtain stringent constraints for individual 
nearby clusters, these rapidly get weaker for more distant 
clusters. In this case, however, one can stack many clusters 
to overcome the loss of signals from each cluster. Reimer 
et al. (2003) took this approach for the EGRET analysis, 
and obtained an improved upper limit to the average flux of 
50 nearby clusters. We argue that the flux is not a physical 
quantity because it depends on distance and therefore distri- 
bution of sources. We should instead convert this improved 
flux limit to constraint on more physical quantities such as 
7-ray luminosity. Here, we examine the GLAST constraints 
on Xp(7?5oo) obtained by stacking clusters from the whole 
sky and in several redshift intervals. As we consider rather 
distant clusters, they are all treated as point sources. 



4.1 Stacking 7-ray signals from galaxy clusters 

^.1.1 Formulation and models 



The number of clusters with M > A/th 
and Z2 is given by 

iVci = / Az— / dM—{M,z), 



between redshifts z\ 



(4) 



where dV^ is the comoving volume element, drih/dM is the 
halo mass function (comoving number density of dark mat- 
ter halos per unit mass range); the former can be computed 
given cosmological parameters, and for the latter we use the 
following parameterization: 



dnh 



= Aj 



Mi; 



■ exp 



f-llncr^' + Bjl'-' 



d In (T 
dMisoii 



(5) 



dMlsOn 

where Pc is the critical density of the present Universe, 
(T(Mi80m, z) is a standard deviation for distribution function 
of linear over density, Aj = 0.315, Bj = 0.61, and ej — 3.8 
(Jenkins et al. 2001). Here we note that i\fi80m is defined as 
an enclosed mass within a given radius, in which the average 
density is 180SlmPc(l + z)^ ■ 
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Figure 2. (a) Cluster mass function as a function of M200 at 
several redshifts. Threshold mass M^^ corresponding to Tt^ = 
3, 5 keV is shown as vertical lines, (b) Cluster mass function 
dnij/d In A/200 multiplied by M|qq((X -Fxp)i in arbitrary units. 
Line types are the same as in (a). 



We give the threshold mass Mii-i{z) in terms of threshold 
temperature Tth based on the observed mass-temperature 
relation: Af200 = 10^^ Mq{T ^ieVf^Eiz)-^ (Voit 
2005). This is because the efficiency of large-scale SZE clus- 
ter surveys relies mainly on cluster temperature regardless 
0/ cluster redshifts. Note that this relation is between tem- 
perature and mass Af200, which is within a radius -R200- Here 
we use the prescription of Hu & Kravtsov (2003) for the con- 
version of different mass definitions, Af200 and Afi80m with 
assumed concentration parameter Cv = 3. For the threshold 
temperature, we adopt Tth = 3 and 5 keV. Fig. 2(a) shows 
the mass function as well as threshold mass corresponding 
to Tth, at various redshifts. In Table 2, we list values of A'^ci 
after integrating equation (4) , for several redshift ranges and 
different Tth. 

The average flux of 7-rays from these clusters is 



^--^ -hiy^t [jM^^iM^z)F.^M,z), (6) 
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Table 2. GLj45T sensitivities to Xp(Rzoo) and yp(-R50o) by stacking N^i clusters above threshold temperature Ttji at given redshift 
ranges, for ap = 2.1, /? = 0, and -Emin = 1 GeV. 
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Figure 3. Relation between 7-ray luminosity (above 100 MeV) 
and cluster mass M20O: for several nearby clusters and for the 
parameters Xp = 0.1, Qp = 2.1, and /3 = 0. Filled (open) points 
are for cooling flow (non-cooling flow) clusters. The solid line is 
the L-y oc M|qq profile that fits the data quite well. 



where Fx^ {M, z) is the 7-ray flux from a cluster of mass M 
at redshift z, given Xp. The flux from each cluster above 



E„ 



is written as 



Fx,{M,z) 



l + z 



dE q^{E,r\M), (7) 



where dL is the luminosity distance, dVd represents the clus- 
ter volume integral, and q.y is the volume emissivity given 
by equation (1) or (Al). 

We then quantify the mass dependence of this flux 
Fxp{M, z). In the case of the isobaric model (/3 = 0) with a 
fixed Xp, the 7-ray luminosity scales as ICM number den- 
sity times energy density, i.e., L^y oc XpUnpth oc XpTV^T. 
On the other hand, luminosity of X-rays due to the ther- 
mal bremsstrahlung process scales as Lx oc n'^T^^^. There- 
fore, there is a relation between 7-ray and X-ray lumi- 
nosities as follows: L^/Lx oc XpT^^^. In addition, there 
are empirical relations between X-ray luminosity and clus- 
ter mass, Lx oc M200, and also between gas temperature 
and mass, T oc mIIIe^/^{z) (Voit 2005). Thus, combin- 



ing these three and assuming that Xp is independent of 
mass, we obtain a scaling relation oc XpM2iME^^^{z). 
In Fig. 3, we show predicted 7-ray luminosity as a function 
of cluster mass (inferred from temperature) , for several well- 
measured nearby clusters (taken from tables in Pfrommer & 
EnBlin 2004) with the parameters Xp = 0.1, Op = 2.1, and 



0^(> 100 MeV) 



7.6 X 10**Xp(M2oo/1015/i-1Mq)21 erg 
s^^ for clusters at z ~ 0, shown as a solid line in Fig. 3. 
When we compute the 7-ray flux Fx^ (-^k) (or equivalently 
luminosity) from clusters with a given mass M, we adopt 
this mass-luminosity relation as a model for average clus- 
ter, and scale E^^^{z) for high-redshift clusters. 

Fig. 2(b) shows the mass function weighed by the mass 
dependence of the flux (in arbitrary unit). This quantity 
represents which mass scale dominates the average flux at 
each redshift. From this flgure, one can see that clusters 
with K'hoQ ~ 3 X 10^'* Mq most effectively radiates 7-rays 
in the low-redshift Universe, but the distribution is rather 
broad for ~1O"-1O^'^M0. If we adopt Tth = 5 keV, then the 
clusters around the threshold mass are the more dominant 
contributors to the average flux. 



4.1.2 GLAST constraints on Xp 

The average flux of the stacked clusters (equation 6) is then 
compared with the corresponding GLAST sensitivity. 



Fii„ 



(8) 



where is the sensitivity to each cluster given as the 

thick dashed line in Fig. 1 (for a point-like source). To de- 
rive constraints on Xp from the stacked image, we solve 
^st,Xp = Jst.iim for Xp. Throughout the following discus- 
sion, we adopt /3 = 0, Op = 2.1 and iJmin = 1 GeV, around 
which the 7-ray yields are maximized compared with the 
point-source sensitivity (Fig. 1). In addition, the pixel num- 
ber with this threshold (47r divided by PSF area; 6 x 10*) 
is large enough to minimize the positional coincidence of 
multiple clusters (compare with A'ci's in Table 2). 

We summarize the results in Table 2. We flnd that the 
limits are as strong as Xp < 0.16 (0.23) for 0.1 < z < 0.15 
(0.15 < z < 0.25). The sensitivities improve for larger Tth, 
because the smaller cluster number is compensated by the 
strong mass dependence of the flux. The constraints on Xp 
degrades rapidly with redshift. Table 2 also shows GLAST 
sensitivities for Yp, which is almost twice as stringent as 
those for Xp in the case of Op = 2.1. We discuss implications 
of this result for Yp in Section 5 in details. 
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The current X-ray catalog covers clusters ai z < 0.2 for 
Tth = 5 keV (Bohringer et al. 2001). The GLAST da.ta. could 
thus immediately be compared with this low-redshift cata- 
log. At higher redshifts, the South Pole Telescope would find 
many clusters with T > 3 keV using SZE; but since it covers 
~10% of the whole sky, the limits would become ~3 times 
weaker than those in Table 2. The Planck satellite, on the 
other hand, would yield all-sky SZE catalog of very massive 
clusters; we find that the limits for Tth ~ 8 keV clusters are 
nearly identical to those for Tth = 5 keV systems. 

In addition to probing its redshift evolution, the stack- 
ing approach is also useful for studying cosmic-ray compo- 
nent in nearby low-mass clusters, and the dependence of Xp 
on cluster mass. Although individual clusters are not bright 
enough, cluster mass function predicts that there are a num- 
ber of such low-mass clusters, which should help improve the 
GLAST sensitivity. 



4.2 Extragalactic 7-ray background 

Another avenue to constrain the universal average of Xp 
is to use the extragalactic 7-ray background (Sreekumar et 
al. 1998), because galaxy clusters would contribute to this 
background intensity to a certain extent. Their contribution 
is quantified as 

'^=1 ^^did^lI/^^dAl^^'^)^-^^^'^)' 

which is quite similar to equation (6). Adopting the same 
models for dwh/dM and Fx^ as in Section 4.1, and using 
Qp = 2.1, /3 = 0, and Emin = 100 MeV, we obtain 

I^{> 100 MeV) = 4 X 10'^ Xp cm"^ s"^ sr"\ (10) 

Even with Xp = 1, this is much smaller than the measure- 
ment by EGRET: 10"^ cm"^ s"^ sr"^ (Sreekumar et al. 
1998). This indicates that cosmic-ray processes in galaxy 
clusters are very unlikely to contribute to the 7-ray back- 
ground flux significantly, especially because it requires a very 
large value for Xp, which is already excluded by EGRET 
for some of nearby clusters. This result is consistent with 
the previous studies such as Colafrancesco & Blasi (1998). 
Hence, we conclude that the stacking method using resolved 
clusters introduced in Section 4.1 would provide much more 
stringent constraint on Xp than the approach using extra- 
galactic 7-ray background. 

However, we here mention a few possibilities that may 
render this approach more viable in the near future. Soon 
after launch, GLAST should start resolving many point 
sources (mainly blazars) that are now contributing to the 
background flux. Furthermore, using angular power spec- 
trum of the 7-ray background map might enable to dis- 
entangle the origin (Ando & Komatsu 2006; Ando et al. 
2007b). In addition, there is a claim that the measured 7- 
ray background flux is dominated by the Galactic foreground 
even at high latitude, and that there is no certain measure- 
ment of truly extragalactic component (Keshet, Waxman, & 
Loeb 2004). In any of the cases above, the contribution from 
galaxy clusters might be found to be significantly smaller 
than the current observed flux, which would be useful to 
constrain Xp at higher redshifts. 
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Figure 4. Relation between ratios of pressure {Yp = Pp/Pth) and 
energy density (Xp = Pp/pth) plotted as a function of spectral 
index Qp of cosmic rays (solid line). Dotted line is the linear fit 
Yp/Xp = 0.5(ap - 1). 



5 X-RAY AND SZE CLUSTER MASS 
ESTIMATES 

Future 7-ray observations of galaxy clusters will have the 
potential to place tight constraints on the nonthermal pres- 
sure provided by cosmic rays. These forthcoming 7-ray con- 
straints will, in turn, provide important handle on system- 
atic uncertainties in the X-ray and SZE cluster mass esti- 
mates based on the hydrostatic equilibrium of the ICM. The 
hydrostatic mass profile of a spherically-symmetric cluster 
is given by 

where M{< r) is the mass enclosed within radius r, pg is 
the gas density, and Pth and Pnt are the thermal and the 
nonthermal contributions to the pressure. The thermal gas, 
measured directly with current X-ray and SZE observations, 
provides a significant fraction of the total pressure support. 
The contribution of the nonthermal pressure, on the other 
hand, is customarily assumed to be small (< 10%) outside of 
a cluster core (see e.g., Nagai, Kravtsov, & Vikhlinin 2007b), 
and it is often ignored in the hydrostatic mass estimates 
based on X-ray and SZE data. The cosmic-ray pressure, if 
present, is a potential source of systematic bias in the hy- 
drostatic mass estimates of clusters (e.g., Ensslin et al. 1997; 
Rasia et al. 2006; Nagai et al. 2007a, and references therein). 

In equation (11), a directly relevant quantity is pres- 
sure gradient rather than energy density Xp that we mainly 
discussed until this point. Currently, it is not possible to in- 
fer both pressure and its radial profile, and here, we simply 
assume that the cosmic-ray pressure profile is the same as 
that of thermal pressure. In this case, one needs to relate 
Xp to Yp. If the cosmic rays are dominated by relativistic 
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component, then equation of state would be Pp = Pp/3. 
On the other hand, for nonrclativistic thermal gas, it is 
Pth = 2pth/3. Thus, we expect Pp/Pth = (l/2)(pp/pth) = 
Xp/2. More precisely, wc can obtain the equation of state for 
cosmic-ray protons by numerically integrating the following 
expressions: 

pp = ^ dp /p(p) (^^p2 _)_^2 _ ^ (12) 

where /p(p) oc p"™" is the differential number density dis- 
tribution. In Fig. 4, we show a correction factor between 
the pressure ratio Yp and Xp, as a function of spectral 
index Qp. This relation is well fitted by a linear formula 
Yp/Xp = 0.5(ap — 1) as shown as a dotted line in Fig. 4; 
the deviation is only ~0.3% at Qp = 2.7. As expected, for 
ap close to 2, the ratio is about 0.5. Therefore, the expected 
sensitivity of GLAST lov Yp would be stronger than that for 
Xp given in Table 1 and as explicitly shown in Table 2. For 
Qp = 2.1, GLAST sensitivities to Yp based on the cluster 
stacking method are 5%, 9%, and 13% at 0.05 < z < 0.10, 
0.10 < z < 0.15, and 0.15 < z < 0.25, respectively. Note, 
however, that the conversion between Yp and Xp depends 
on ap, for which lACT measurements would be essential. 

Observational constraints on Xp = (pp) / (pth) is also 
sensitive to any non-negligible small-scale structure in the 
ICM. When gas clumps, it has density higher than the local 
average, (pth). If it is not resolved and masked out, the local 
inhoniogcncity in the ICM boosts 7-ray surface brightness 
by a factor of C-, = (ppPth) / (pp) (pth) and X-ray surface 
brightness by Cx = (pth) / (Pth)^, while leaving SZE signal 
(which is linearly proportional to pth) unaffected by dumpi- 
ness. A joint 7-ray-|-X-ray constraints on Xp based on a 
smooth model is generally biased by a factor C-^/Cx, which 
could be greater or less than 1 depending on the relative 
size of and Cx-^ A joint 7-ray -l-SZE constraint on Xp, 
on the other hand, is biased high by a factor C-y. Recent 
cosmological simulations of clusters that include cosmic-ray 
physics indicate jagged shape for the Xp(t) profile, which 
implies a large clumping (Pfrommer et al. 2007). These 
simulations are potentially useful for estimating the values 
of C-y, which would be important for interpretation of Xp 
in case of detection of cluster signals with upcoming 7-ray 
experiments. In absence of these constraints, observational 
constraints on Xp should be taken as an upper limit. 

Recently, Mahdavi et al. (2007) performed a compari- 
son between masses estimated with weak gravitational lens- 
ing and using the assumption of hydrostatic equilibrium, 
and showed that the latter masses are typically biased to 
be lower by 20%. This result might indicate presence of the 
nonthermal pressure component. Upcoming 7-ray measure- 
ments of galaxy clusters could thus provide useful informal 
tion on the origin of this mass discrepancy. Turbulence and 

^ Current X-ray observations with superb spatial resolution and 
sensitivity are capable of detecting the prominent clumps that 
contribute significantly to the X-ray surface brightness. A com- 
parison of recent X-ray and SZE measurements indicate that the 
X-ray clumping factor is very close to unity (1 < Cx ^ 1-1) in 
practice (LaRoque et al. 2006). 



magnetic fields are also potential sources of bias in X-ray 
and SZE cluster mass estimates. Recent numerical simula- 
tions of cluster formation indicate that sub-sonic motions of 
gas provide nonthermal pressure in clusters by about ~10% 
even in relaxed clusters (e.g., Rasia et al. 2006; Nagai et al. 
2007a, and references therein). Most cluster atmospheres are 
also magnetized with typical field strengths of order a few 
ptG out to Mpc radii (Carilh & Taylor 2002; Govoni & Fer- 
etti 2004), but this would only give negligible contribution 
to the total pressure support. 



6 INVERSE-COMPTON SCATTERING FROM 
NONTHERMAL ELECTRONS 

6.1 Secondary electrons from pion decays 

Until this point, we have neglected the contribution to 7- 
rays from relativistic electrons and positrons produced from 
decays of charged pions. Those charged pions are produced 
by the proton-proton collisions just as tt^'s that decay into 
7-rays. Thus, as long as the cosmic-ray protons exist, there 
should also be relativistic e* component associated with 
them. GeV 7-rays would be produced by IC scattering of 
CMB photons due to such a 'secondary' leptonic compo- 
nent. In this subsection, wc show the expected IC flux to 
compare it with the flux from vr'^ decays, and argue that the 
former is indeed negligible, justifying our earlier treatment. 

Unlike protons, leptons can cool quickly by synchrotron 
radiation and IC scattering. Energy distribution of these 
electrons (positrons) after cooling is obtained as a steady- 
state solution of the transport equation, which is 

n,{E,,r)= / dS:Qe(£;:,r), (14) 

where Qc is the source function of injected electrons. For the 
energy-loss rate Eg, the dominant interaction would be syn- 
chrotron radiation and IC scattering of CMB photons, i.e., 
—Ee oc {U^ + Ucms)E'^, where J7b and f/cMB are the energy 
densities of magnetic fields and CMB. If the injection spec- 
trum is power law, Qc oc E'"" , then equation (14) states 
that the spectrum after cooling would be Ue oc E'""'^, 
steeper by one power. 

Once we know the electron distribution wc can unam- 
biguously compute the IC spectrum after scattering CMB 
photons. In addition, in the case of the secondary electrons, 
we can compute the source Qe relatively well given the spec- 
trum of cosmic-ray protons. In Appendix B, we summarize 
fitting formula that wc use, given by Dolag & Enfilin (2000) 
and Pfrommer & Enfilin (2004). Looking at equation (14), 
in order to get the electron distribution after cooling, we 
also need to know magnetic field strength B in the clusters 
that is relevant for synchrotron cooling. The estimates of 
B range '^O.l 10 pG (Clarke, Kronbcrg, & Bohringcr 2001; 
Fusco-Femiano et al. 2004; Rephaeli, Gruber, & Arieli 2006), 
while the CMB energy density corresponds to equivalent 
field strength of Bomb = 3.24(1 -I- z)'^ pG. Thus, unless B is 
larger than or comparable to Bcmb everywhere in the clus- 
ter, the synchrotron cooling would not be significant, as the 
energy loss is proportional to B^ + Bomb- We here assume 
S = to obtain the maximally allowed IC flux. 

In Fig. 5(a), we show flux of IC 7-rays from secondary 
leptons, compared with direct 7-ray flux from 7r° decays. 
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Figure 5. (a) Flux of 7-rays from tt^ decays with Xp = 0.1 
(dotted), IC scattering due to secondary electrons (dashed) as 
a function of minimum energy Bminj for Qp = 2.1, /3 = 0, and 
B = 0; total flux is indicated as a solid curve, (b) Fractional 
contribution of IC scattering to the total 7-ray flux, i^ic/^tot> 
for a = 2.1 (solid), 2.4 (dotted), and 2.7 (dashed). 



assuming Xp — 0.1, Qp = 2.1, and (3 — 0. Fig. 5(b) shows 
the fractional contribution of the IC processes for various 
values of Qp. These figures show that even in the case of very 
week magnetic fields to reduce the electron energy losses, the 
IC processes give only sub-dominant flux in the GeV energy 
range relevant for GLAST. The fractional contribution of 
the IC emission to the total 7-ray flux, which is independent 
of Xp, is smaUer than 20% for E^in = 100 MeV and a = 
2.1. For a steeper proton spectrum (a > 2.1), the fractional 
contribution become considerably smaller. Bremsstrahlung 
process due to the same electrons and positrons is even more 
suppressed (Blasi 2001). We thus conclude that the IC and 
bremsstrahlung 7-ray emission by secondary electrons are 
sub-dominant for the realistic range of parameters. 



6.2 Primary electrons by shock acceleration 

Whenever the shocks are generated, both ions and electrons 
are accelerated. Thus, one expects that the IC scattering off 
the CMB photons due to such primary electrons would also 
contribute to the GeV-TeV 7-ray flux to a certain extent 
(Loeb & Waxman 2000; Totani & Kitayama 2000; Waxman 
& Loeb 2000; Gabici & Blasi 2004). If this process dominates 
the 7r° decays in 7-ray energy band, then the constraints on 
Xp will be directly affected in case of detection. However, 
there are difficulties for this mechanism to work efficiently 
in many clusters. 

As electrons lose their energies via radiation much more 
rapidly than protons, clusters would be bright with this 
mechanism during only a limited period after injection. For 
example, the radiative cooling time scale for 10 GeV elec- 



trons is ~10* years, which is much shorter than typical clus- 
ter age. By the same reason and also comparing the spatial 
intensity distribution, it is unlikely that synchrotron radia- 
tion from these primary electrons is responsible for the ob- 
served radio halo emissions (e.g., Blasi, Gabici, & Brunetti 
2007). 

It might still be possible to overcome these difficul- 
ties if these electrons are continuously reaccelerated in situ 
through the second order Fermi mechanism (Schlickeiser, 
Sievers, & Thiemann 1987; Tribble 1993; Brunetti et al. 
2001; Petrosian 2001). In this case, however, the spectrum 
of electrons has typically a cutoff at the Lorentz factor of 
< 10^ . This property, while explains spectrum of radio halo 
of Coma quite well (e.g., Reimer et al. 2004), would restrict 
the 7-ray flux in the GeV region due to the IC scattering 
and bremsstrahlung. In Fig. 1, we show the upper bound on 
these components in the case of Coma cluster as a dotted 
curve, taken from Reimer et al. (2004). 

In consequence, as long as Xp is more than a few 
per cent, it would be unlikely that the primary electrons, 
whether they are directly injected or continuously reacceler- 
ated, dominate the GeV 7-ray flux, at least in a large frac- 
tion of clusters. Even though primary electron component 
dominated the detected flux, the shape of 7-ray spectrum 
would be very different from Tr'^-decay component especially 
at low energies, which could be used as a diagnosis tool; 
this difference comes from the kinematics of tt" decays. The 
GLAST energy band ranges down to ~20 MeV, which is 
especially important characteristic for that purpose. More- 
over, observations in lower frequency bands such as radio, 
EUV, and hard X-rays, are also important, because these 
emissions are understood as synchrotron radiation (for ra- 
dio) and IC scattering (for EUV and hard X-rays) from non- 
thermal electrons. 



6.3 Secondary leptons from ultra-high energy 
cosmic-ray interactions 

If protons are accelerated up to ultra-high energies such as 
> lO'^* eV in galaxy clusters, which may be plausible, these 
protons are able to produce e* pairs through the Bethe- 
Heitler process with CMB photons: P7cmb — > pe"'"e~. These 
high-energy e^ pairs then IC scatter the CMB photons, pro- 
ducing GeV-TeV 7-rays (Aharonian 2002; Rordorf, Grasso, 
& Dolag 2004; Inoue, Aharonian, & Sugiyama 2005). In this 
case, the IC photons might dominate the vr^-decay 7-rays by 
many orders. 

However, this mechanism is extremely uncertain, de- 
pending heavily on the maximal acceleration energy of the 
protons. This is especially because the threshold energy of 
the Bethe-Heitler process is ~10^^-10^® eV, and it IS un- 
clear whether the magnetic fields are strong enough to con- 
fine these ultra-high energy protons for cluster ages. Even if 
the detected 7-rays are dominated by this mechanism, the 
spectrum would be quite different from the vr^-decay 7-rays 
and should be easily distinguishable (e.g., Inoue et al. 2005). 
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7 CONCLUSIONS 

We investigated the capability of tlie current and future 7- 
ray detectors such as GLAST and lACTs for constraining 
the cosmic-ray pressure contribution to the ICM. 
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(i) We showed that the upcoming GLAST measurements 
can be used to place stringent upper limits, 0.5-5%, on the 
ratio of energy densities of the cosmic rays and thermal gas, 
Xp, for several nearby rich clusters. These limits are fairly 
insensitive to the assumed energy spectrum or the radial dis- 
tribution of the cosmic-ray protons for a reasonable range 
of models. We showed that lACT sensitivity to Xp is not as 
stringent as GLAST, but lACTs provide useful constraint 
on spectral index Op, which in turn provide important con- 
straints on the acceleration mechanisms of cosmic rays. 

(ii) The stacking method offers a powerful technique to 
probe the cosmological evolution of Xp and Fp with upcom- 
ing 7-ray observations. Using the latest cosmological models 
such as halo mass function and phenomenological relations 
that reproduce observed cluster properties, we showed that 
one-year all-sky survey with GLAST can place tight limits 
(Fp < 10%) on the evolution of mean cosmic-ray pressure 
in clusters out to fairly high redshift [z < 0.25) by stacking 
signals from a large sample of known clusters. These con- 
straints will correspond to an upper limit on the systematic 
uncertainties in the X-ray and SZE cluster mass estimates, 
due to nonthermal pressure provided by cosmic rays. In ad- 
dition, since the halo merger rate is expected to increase 
with redshift (e.g., Gottlober, Klypin, & Kravtsov 2001) and 
such mergers can boost 7-ray signals (Pfrommer et al. 2007), 
the technique may provide insights into the relation between 
cosmic-ray energy density and merger activities. The same 
approach will also enable one to probe cosmic-ray popula- 
tions in low-mass clusters. 

(iii) We also evaluated the cluster contribution to the ex- 
tragalactic 7-ray background using the latest models, and 
showed that even with Xp = 1, the contribution is only 
about 4% of the measured flux. This indicates that this ap- 
proach would not currently be very helpful to constrain Xp, 
but might become more useful in the future if a significant 
fraction of the background flux were resolved. 

(iv) We showed that 7-rays due to IC scattering by both 
the primary and secondary electrons are likely sub-dominant 
relative to the 7-rays from 7r° decays in most of the clusters. 
We flnd that the fractional contribution of the IC flux by 
secondary electrons never exceeds ^20% for a reasonable 
range of parameters, independently of Xp. The contribution 
from the primary electrons will also be suppressed in many 
clusters, because either they cool very fast after injection or 
they cannot be accelerated up to very high energies in the 
reacceleration models. Moreover, multi-wavelength observa- 
tions in radio, EUV, and hard X-ray wavebands will provide 
independent constraints on nonthermal electrons in clusters 
(e.g., Reimer et al. 2004), and such a consideration shows 
that the expected 7-ray flux from the primary electrons is 
indeed sub-dominant as long as Xp > 0.02 (Fig. 1). Even if 
these components were dominant in some clusters, the shape 
of 7-ray spectrum should provide diagnostics of the origin. 
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APPENDIX A: 7-RAY EMISSIVITY FROM n° 
DECAYS 

Equation (1) has a very clear structure including several rele- 
vant physics, ranging from cosmic-ray distribution np{Ep, r) 
to 7r°-production cross section dapp/dE^- This integral is no 
difficult, and indeed, Pfrommer & Enfilin (2004) gave a sim- 
ple fitting form for that as follows: 



q-f{E,r) = fTeff criH (r)^ 



GeV 3a-y \ GeV 



2E 



(Al) 



where — Op is the asymptotic spectral index of 7-rays 
that is the same as that for protons, 5^ = 0.14q:^^ '' + 0.44, 
^ = 2 is a constant pion multiplicity, and 



cTeff = 32(0.96 + e*-4-2 -*«7 ) mb 



(A2) 



is the effective inelastic pp cross section. This reproduces 
results of numerical computations of hadronic processes as 
well as accelerator data quite well. 

As 7r°'s are produced by collisions between nonthermal 
cosmic-ray ions and thermal ICM nucleons, 7-ray emissiv- 
ity is proportional to the product of ICM density nu and 
number density of cosmic rays. The latter quantity is ef- 
fectively characterized by hp and this is given by requiring 
that the fraction Xp of the thermal energy density pth goes 
to cosmic-ray energy density: 



Xp{r)pth{r) 



hp{r)mpC"' /w-pC^V 
2(ap - 1) I, GeV J 



(A3) 



where B is the beta function, appearing when we integrate 
kinetic energy of each proton weighed by the momentum 
distribution function, and 



Pth(r) = - 1 + 



1 - 3XHe/4 
1 - XHe/2 



ne(r)feBTe(r), 



(A4) 



with fee the Boltzmann constant, Xue = 0.24 is the primor- 
dial mass fraction of *He, and the electron density rie and 
temperature are well measured with X-rays. 



APPENDIX B: INVERSE-COMPTON 
SCATTERING FROM SECONDARY 
ELECTRONS 

Hadronic collisions also produce charged pions that eventu- 
ally decay into electrons and positrons. These leptons, hav- 
ing relativistic energies, can up-scatter the CMB photons 
into GeV energies. Since the physics of IC scattering is well 
established (Rybicki & Lightman 1979) and pion produc- 
tion due to pp collisions are measured in laboratories, this 
process can be described with relatively small ambiguity. 
Electron distribution function after radiative cooling is 



ne{Ee,r) = 



hc{r) f E, 



GeV VGcV 



he{r) 



2^7rl6-'"'=-^' o-effmeC* nii{r)hp{r) 



(Bl) 



(B2) 



Oe - 2 aT GeV B(r)2 + Bg^g ' 

where Qc = Qp + 1, itt is the Thomson cross section, and 
BcMB = 3.24{1+ z)^ fxG. Emissivity of IC scattered photons 
is given as 



qic{E,r) = q{r)fic{ae) 



f meC 



,2 \ l-ae 



(B3) 



© 2008 RAS, MNRAS 000, 1-12 



12 S. Ando and D. Nagai 



q{r) = 
/ic(ae) = 



87i-Vgne(r)(fcBrcMB)'' 
2"''+3(a2+4Qe + ll) 



(ae + 3)2(a 

^ /Oe + 5 

X r 



5)(qc + 1) 
^ , ae + 5 



(B4) 



(B5) 



where a„ = (qo — l)/2, Tc is the classical electron radius, 
TcMB = 2.7 K is the CMB temperature, F is the F-function, 
and C is the Riemann (^-function. For more detailed discus- 
sions, see Dolag & EnfiUn (2000) and Pfrommer & EnfiUn 
(2004). 



APPENDIX C: GLAST SENSITIVITY FOR 
EXTENDED SOURCES 

Flux sensitivity of GLAST-LAT to a point-like source is 
shown as a thick dashed line in Fig. 1. The sensitivity for 
an extended source is different; in this section, we derive it 
using a simple argument. 

The dominant background is extragalactic 7-ray flux, 
and its intensity depends on photon energy as Ihg oc E~^'^ 
(Sreekumar et al. 1998). At low energies where background 
photon count are larger than 1 (A'bg > 1), the flux sensitivity 
is determined by a criterion such as Mim > a^jN-og, where 
some number a sets significance of detection; hereafter use 
Q = 5. On the other hand, at higher energies where A'bg < 1, 
then the detection simply relies upon photon count from a 
source. 

Thus, for a source like galaxy clusters, there are up to 
four different energy regimes depending on the source ex- 
tcntion. (i) At lowest energies where PSF (or pixel size) is 
larger than the source size (i.e., fJpix > n), the source can 
be regarded as a point-like object. The other three regimes 
are for more energetic photons where source are extended 
(il > fipix); they are where background photon counts are 
(ii) larger than 1 from one pixel (A^bg.pix > 1); (iii) smaller 
than 1 from one pixel but larger than 1 from the entire source 
region (iVbg.pix < Ij.'Vbg > 1); (iv) smaller than 1 from the 
entire source region (iVbg < 1)- For the lowest energy region 
(i), the sensitivity is the same as that for point sources, and 
this corresponds to the regime below ~300 MeV in Fig. 1. 

For the region (ii), the point-source flux sensitivity 
^iim,pix are obtained by a criterion Mim,pix = Sy'iVbg.pix- 
These photon numbers are related to the flux and back- 
ground intensity through Afiini,pix = -Fiim.pix^r and 
A'bg.pix = /bg^TTipix, where AT (effective area times expo- 
sure time) is the detector exposure. The point-source sensi- 
tivity is thus obtained by 



= 5 



/bgQpi 

AT 



(CI) 



A similar argument cam be applied for the flux sensitivity 
for an extended source and we obtain 



Flin, = 5 



AT 



(C2) 



Thus the sensitivity becomes weaker by a factor of 9/50psf, 
compared to that for a point-like source. This is the case for 
the region between 300 MeV and 2 GeV in Fig. 1. 



When the photon energy becomes higher, the back- 
ground count gets smaller. We then consider the region (iii). 
In this case, to obtain the point-source sensitivity, we use 
a criterion of flve-photon detection: Fiim.pix — 5 /AT. As 
we have more background photons (than 1) from the en- 
tire source region, the extended-source flux sensitivity is the 
same as the flrst equality of equation (C2). Combining these 
two, we obtain 



^lim — \/5Flim,pix/bgf2, 



(C3) 



and this is for the region between 2 GeV and 30 GeV in 
Fig. 1. At the highest regime (iv), the region above 30 GeV in 
Fig. 1, the cluster detection is totally relies on photon count 
and independent of background. Therefore, the sensitivity 
for an extended source is the same as that for a point-like 
source. 
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